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INTRODUCTION 

Among three basic states of matter (solid, liquid, gas), 
liquids are least understood from the theoretical point 
of view. The inter-atomic, or inter-molecular, interac- 
tions in a liquid are strong, and therefore strongly affect 
the liquid energy. At the same time, the interactions are 
system-specific, hence the calculation of liquid energy re- 
quires the explicit knowledge of the interactions. For 
this reason, it is argued 1] that no general expressions 
for liquid energy can be obtained, in contrast to gases 
and solids. 

According to current theoretical understanding, liq- 
uids are strikingly different from both gases and solids. 
Indeed, small atomic displacements in a solid make it 
possible to expand the energy in terms of phonons and 
obtain general expressions for solid energy. On the other 
hand, this has largely considered to be impossible to do in 
liquids where atomic displacements are large. Similarly, 
small interactions in gases make it possible to treat inter- 
actions as a perturbation and obtain general expressions 
for the corrections to the system energy. On the other 
hand, this has not been useful for liquids with strong in- 
teractions and solid-like densities. An apt summary of 
this state of affairs, attributed to Landau, is that liq- 
uids "have no small parameter" . Perhaps for this reason, 
liquid heat capacity is not, or is barely, mentioned in sta- 
tistical physics textbooks as well as books dedicated to 
liquids [l-Q. This observation is shared by Granato 0, 
who further comments on the challenge faced by teachers 
to discuss liquid heat capacity in class. 

The experimental behavior of liquid heat capacity is in- 
teresting. The heat capacity per atom decreases from ap- 
proximately Sfce per atom at the melting point to about 
2fcB at high temperature, as witnessed by the data of 
many simple liquids 1^, 0| . Similar behavior is also seen 
in complex liquids [lO|. The decrease of heat capacity 
was observed in molecular dynamics simulations [ll|, [l^ . 
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Liquid heat capacity has also been discussed by consid- 
ering liquid potential energy landscape with inter-valley 
motions fo^. 

Theoretically, liquids have been viewed to occupy an 
intermediate state between gases and solids. Liquids are 
fluid, and therefore might intuitively appear closer to 
gases in terms of their properties. On the other hand, 
unless they are close to critical point, liquids have solid- 
like densities and can support shear waves at high fre- 
quencies. The question of how best to treat real dense 
liquids, i.e. approach them from the gas or solid state, 
has a long history 

Historically, liquids have predominantly been ap- 
proached from the gas state, by developing schemes to 
calculate the interaction energy in addition to gas kinetic 
energy [ll-@|- This approach relies on the knowledge of 
interatomic interactions as well as correlation functions 
13l Il4| . These are not generally available, apart from 
simple model systems, as discussed below in more de- 
tail. When they are available, it is not apparent how 
this approach explains the experimental decrease of heat 
capacity from S/cb to 2k^ at high temperature. 

An alternative approach to liquids was pioneered by 
Frenkel Q, and is based on liquid relaxation time r. r 
is the average time between two consecutive local struc- 
tural rearrangements in a liquid at one point in space. In 
this picture, a liquid is approached from the solid state 
because locally, liquid structure remains unchanged, i.e. 
the same as in a solid, during time shorter than r. An im- 
portant advantage of this approach is that strong interac- 
tions are included in the consideration from the outset. 
This is in contrast to the approach from the gas state 
that attempts to calculate the interaction energy from 
correlation functions and interatomic interactions. 

Frenkel noted that as t continuously increases (if crys- 
tallization is avoided) beyond the experimental time 
scale, a liquid becomes a solid for practical purposes. 
Hence, a solid is different from a liquid only quantita- 
tively but not qualitatively. Frenkel subsequently stated: 
"...the classification of condensed bodies into solids and 
liquids... has a relative meaning convenient for practi- 
cal purposes but devoid of scientific value" [3 . Novel for 
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that time, this view might perhaps come across as some- 
what unusual even at present. We suggest that a possi- 
ble reason for this is that solid-like properties of liquids 
have not become apparent in the traditional approach to 
liquids from the gas state. On the other hand, we pro- 
pose here that as far as their energy and heat capacity 
are concerned, liquids can be understood on the basis of 
their solid-like properties, consistent with Frenkel's gen- 
eral idea. 

We have recently proposed how liquid energy can be 
calculated in Frenkel's approach, by relating liquid en- 
ergy to T [31. We have found good agreement of liq- 
uid heat capacity between theoretical and experimental 
data for mercury. Two important questions subsequently 
arise. First, how well does the theory work for a larger 
number of liquids? Second, the proposed theory was har- 
monic. On other hand, anharmonic effects are known 
to be particularly large in liquids, where the coefficients 
of thermal expansion considerably exceed those in solids 
[l5| . Hence, it is important to extend the theory to the 
anharmonic case and compare it with experiments. 

In this paper, we propose how to include the anhar- 
monic effects in the approach to liquids from the solid 
state. We subsequently calculate heat capacity from both 
harmonic and anharmonic theory, and compare the re- 
sults with experimental data for 5 commonly discussed 
liquids. We find good agreement between theoretical pre- 
dictions and experimental data with no free fitting pa- 
rameters. Finally, we discuss and compare the proposed 
theory to the previous approaches to liquids. 



HARMONIC THEORY 

We begin our discussion with the work of Frenkel Q , 
who provided a microscopic description of Maxwell phe- 
nomcnological viscoelastic theory of liquid flow [l^, by 
introducing liquid relaxation time t. As mentioned ear- 
lier, r is the average time between two consecutive atomic 
jumps in a liquid at one point in space. Each jump can 
approximately be viewed as a jump of an atom from its 
neighboring cage into a new equilibrium position, with 
subsequent cage relaxation. These atomic jumps, or lo- 
cal relaxation events (LREs), give a liquid its ability to 
flow, r is a fundamental flow property of a liquid, and is 
directly related to liquid viscosity 77 as 77 = GqcT (il. [16|. 
where Goo is the instantaneous shear modulus. On tem- 
perature decrease, r increases by many orders of magni- 
tude, reaching 10^ — 10^ s at which point, by convention, 
a liquid forms a solid glass because LREs stop operating 
on a typical experimental time scale. 

In Frenkel's theory, motion of an atom in a liquid con- 
sists of two types: vibrational motion around an equilib- 
rium position as in a solid, with Debye vibration period 
of about td =0.1 ps and diffusional motion between two 
neighboring positions during time t. Therefore, if the 



observation time is smaller than r, the local structure of 
a liquid does not change, and is the same as that of a 
solid glass. In this picture, Frenkel realized that a liq- 
uid should maintain solid-like shear waves, similarly to 
those existing in a solid, at all frequencies lj larger than 
- 0. This prediction was later confirmed experimen- 
tally 17h21||. Longitudinal waves, associated with den- 



sity fluctuations, are considered to be unaffected in this 
picture, apart from different dissipation laws for cj < - 
andtJ>i[2|. 

Basing on the ability to support high-frequency shear 
waves, we have calculated liquid energy as follows [lij . 
In Frenkel's picture, liquid energy is 



E — E,, 



E, 



(1) 



where E,, and Ed are the energies of vibration and diffu- 
sion, respectively. 

The vibration energy consists of the energy of one 
longitudinal mode and two shear modes with frequency 
a; > i. Then, E^ can be written as E^ ^ Ki + Pi + 
Ks{uj > ^) + Ps{uj > where K and P correspond to 
kinetic and potential terms, respectively. Here, the ab- 
sence of shear modes with frequency ut < ^ is expressed 
as Ks{uj < i) = Ps{uj < i) =0, implying the absence of 
restoring forces for low-frequency vibrations. Similarly, 
Ed can be written as Ed — Kd + Pd, where Kd and Pd 
are respective kinetic and potential terms, giving 

E = Ki + Pi+ K,{lo >-) + P^[u>-) + Kd + Pd (2) 

T T 

We now approach a liquid from the solid state (glass) 
where diffusion is absent so that Ed = and E = Ey. At 
certain temperature, shear waves with frequency uj < 
disappear. The oscillatory motion at frequency smaller 
than < i is substituted by diffusional motion tak- 
ing place during time r. The latter motion can be 
thought of as a "slipping" motion at low frequency be- 
cause the restoring force for low-frequency vibrations be- 
comes small (see Figure 1). We note that forces experi- 
enced by diffusing atoms are of the same nature as those 
experienced by vibrating atoms, and are defined by in- 
teratomic interactions that give rise to both Pd and Ps- 
Therefore, the relative smallness of Ps at small frequency 
implies the smallness of Pd as compared to ^^(0; > i) 
(as well as to Pi). In other words, if the interaction en- 
ergy of diffusing atoms, Pd, were large and comparable 
to Ps{uj > i), this would imply strong restoring forces 
and consequently the existence of low-frequency vibra- 
tions. Therefore, Pd can be neglected in Eq. This 
is the only approximation in the theory 15| . We further 
note that the energy of diffusion becomes unimportant at 
low temperature [22|. Indeed, because the time that an 
atom spends in the transitory diffusing state is approx- 

N ' 



imately td, the relative number of diffusing atoms. 
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FIG. 1; As temperature increases, the low-frequency oscilla- 
tory shear motion of the highlighted atom between the two 
dashed lines is lost, as the restoring force for low-frequency 
shear vibrations becomes weak. Instead, the atom "slips" 
and diffuses to another position defined by the solid arrow. 
Schematic illustration. 



is equal to the jump probability, where Nd and N 
are the number of diffusing atoms and the total number 
of atoms, respectively. When r significantly exceeds td, 
the number of diffusing atoms and, therefore, diffusing 
energy, become small and can be ignored. 

We now note that the sum of all kinetic terms in Eq. 
(HI) gives the total kinetic energy of the liquid, K. In- 
deed, K = regardless of how the motion partitions 
into vibrational and diffusional motion. Then, Eq. ([2]) 
becomes 



E = K + Pi+Ps{uj> 



(3) 



It is convenient to re- write Eq. ([3]) using the equiparti- 



tion theorem that implies Pi 



Ml 

2 



E 



2 , where we noted that liquid 
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and K = Ki+K, 

kinetic energy is the same as in the solid (K — ^ j, 
and can therefore be written as a sum of kinetic terms 
related to longitudinal and shear waves. Then, E in 

Eq. © becomes E ^ Ei + ^ + For sub- 

sequent calculations, it is convenient to further write 
Es = Es{uj < ^) + Es{uj > where the two terms 
refer to their solid-state values. Then, E becomes finally 



E^Ei+E^iu:>'-) + ^^^^^^^ (4) 

T Z 



Z2 = (27r?i)-^' / exp(-^f^(p; 



^Wt) dpdq (5) 



2 ' uj^qj) j dpdq 



where = 7 , '-^h and ujgi are frequencies of longitudinal 
and shear waves, N is the number of atoms and N' is the 
number of phonon states that include longitudinal waves 
and transverse waves with frequency uj > Here and 
below, /cb = 1. 

Integrating, we find 



N 



N 



T 



1=1 



/ 2N 

n 



(6) 



where A'^i is the number of transverse modes with uj > loq. 

In the harmonic approximation, frequencies uju and wu 
are considered to be temperature-independent, in con- 
trast to anharmonic case discussed in the next section. 
Then, Eq. ([5]) gives the liquid energy in harmonic ap- 
proximation E ^ T'^-§p\nZ = NT + NiT. Ni can be 
calculated using the quadratic density of states in the 
Debye model, as is done in solids IJ. The density of 
states of shear modes is 34(0;) = -^^w^, where uj„it is 
Debye frequency of shear modes (wmt ~ wd) and we 
have taken into account that the number of shear waves 

is 2N. Then, A^i = / gt{uj)duj = 27V f 1 - f ^ ' 

To calculate the last term in Eq. (|4|), we note that 
similarly to Et{uj > 1/t) = NiT, Et{uj < 1/r) = N2T, 
where N2 is the number of shear modes with uj < luq. 

Because N2 = 2N - N , N,, = 2N { \ . Then, E 



Ni + ^)r, giving finally 



E = NT\'S 




(7) 



According to Eq. (I7|), the liquid energy is 3NT as in 
a solid when r ;3> td. This gives the heat capacity per 
atom, Cv = j^^j Cy = 3, consistent with experimental 
results d, II]. When r — td at high temperaturej__Eq. 
([7]) predicts c„ = 2, as in the experimental data [T 



The first two terms can be calculated in the same way 
as is done in the harmonic theory of solids [1], except 
we separate longitudinal and shear waves in the partition 
function and account for shear waves with frequency uj > 
- only. Let Z2 be associated with the first two terms in 
Eq. (gl). Then, Z2 is: 



ANHARMONIC THEORY 

In the harmonic approximation, temperature depen- 
dence of frequencies in Eq. ([6]) is ignored. This may be 
a good approximation for phonons in solids at either low 
temperature or in systems with small thermal expansion 
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coefficient a. In liquids, on the other hand, anharmonic- 
ity and associated thermal expansion are large (sij . and 
need to be taken into account. In this case, applying 
E = T^-^lnZ to Eq. © gives 



the number of shear modes with lo < luq defined in the 
previous section. Adding this term to Eq. pO)) and using 
Ni and N2 calculated in the previous section gives the 
anharmonic liquid energy: 



E2 = NT-T^y —^+N,T-T^y 

^ oJii dT ^ ujsi dT 



(8) 



which corresponds to the first two terms in Eq. 

The anharmonic effects lead to the decrease of fre- 
quencies with temperature, resulting in system soften- 
ing [sH. We therefore need to calculate and 
in Eq. ([8]). These can be calculated in the Griineisen 
approximation, by introducing the Griineisen parameter 
7 = — (^^) , where are frequencies in harmonic 
approximation [22|. 7 features in the phonon pressure 
Pph = — (^)^, where F is free energy in the harmonic 
approximation. F = —TlnZ can be calculated from Eq. 
([6]), giving 




(11) 



At low temperature when r exceeds td, Eq. (jlip gives 
E = 3N{T + and c„ is 



Cv 



3(1 + aT) 



(12) 



Eq. (jl2p explains why experimental Cy of liquids ex- 
ceed the Dulong-Petit value just above the melting point 
In the next section, we compare the calculated and 
experimental Cv in the entire temperature range. 
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Fo 



T 



i=l 



(9) 



Calculating Pph — — (f^)^ and introducing the above 
Griineisen parameter gives Pph — -y-{N + Ni). Then, 
the bulk modulus due to the (negative) phonon pres- 
sure is B„h - - w ^ ATA „:.Ar.„ ( ^ 
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( dBph\ 



3p/l = V gy = -tr(A^ + Ni), glVmg y-Q^ 

y\j.^ I A'^i). We now use the macroscopic definition 
of 7, 7 = ^"-B — 



33 • Here, B — Bq + Bph is the total 



bulk modulus, Bq is zero-temperature bulk modulus, a 
is the coefficient of thermal expansion and Cy is constant- 
volume heat capacity. For C^,, we use its harmonic value, 
Cy = N + Ni from Eq. dU, because ^ and ^ in 
Eq. ([8|) already enter as quadratic anharmonic correc- 
tions. Then, -a{BQ + Bph). For smaU aT, 

-T, 



as is often the case in experiments, this implies B oc 
consistent with the experimental data ^32| . 

We note that experimentally, B linearly decreases with 
T at both constant volume and constant pressure [3^ . 
The decrease of B with T at constant volume is due to 
the intrinsic anharmonicity related to the softening of in- 
teratomic potential at large vibrational amplitudes; the 
decrease of B at constant pressure has an additional con- 
tribution from thermal expansion. 

Assuming uP' a Bq -\- Bph and combining it with 
'dB^\ ^ + from above gives ^ 

Putting it in Eq. ([8]) gives: 




COMPARISON WITH EXPERIMENTAL DATA 

We now compare the predictions of Eqs. JT]) and pT|) 
with the experimental Cy. Simple liquids [a @| offer a 
good test case for a number of reasons. First, they are 
frequently discussed in the literature. Consequently, the 
data of both Cy and viscosity are available for simple 
liquids, partially because their melting temperatures are 
fairly low, simplifying the measurements. Second, these 
liquids become low- viscous close to the melting point, in 
contrast to, for example, liquid B2O3, Si02 and so on. 
Consequently, the decrease of Cy, predicted by Eqs. ^ 
and (llip to operate around t « td, can be observed in 
a convenient and accessible temperature range. Third, 
it is important to see that our theory works in a wide 
temperature range. We have taken experimental Cy for 
liquid Hg, In, Rb, Cs and Sn where Cy has been measured 
in a wide range of temperatures, from 200 to 1200 K. 

To the best of our knowledge, apart from liquid Hg, 
this is the first attempt to analytically calculate Cy of the 
above liquids and compare them with experiments. 

To calculate Cy from the theory, we write r = ^ in 
terms of liquid viscosity using the Maxwell relationship 
T = Then, harmonic and anharmonic Cy in Eq. ([7]) 
and Eq. (jlip are given by the two equations below: 




(13) 



E2 ^{N + Ni){T 



(10) 



2 

in the same way, giving (T + §T^), where N2 is 



The last term in Eq. (jH) , ^^'^'^^^/'^i ^ ^an be calculated 




T + 




(14) 



We have taken viscosity data from Refs [25|-|27| and 
fitted them to the form of the Vogel-Fulcher-Tammann 
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FIG. 2: (Color online). Experimental and theoretical for 
liquid Hg. Experimental is shown in square symbols. Cal- 
culated harmonic and anharmonic c„ are shown in the bottom 
and thick black curve, respectively. tdGoo ~ 5.5 ■ 10^* Pa-s 
and 6.272-10"'* Pa-s for theoretical harmonic and anharmonic 
c„, respectively, a — 1.60 - 10"'* K~*. 



FIG. 3: (Color online). Experimental and theoretical for 
liquid In. Experimental Ci, is shown in square symbols. Cal- 
culated harmonic and anharmonic Cv are shown in the bottom 
and thick black curve, respectively. tdGoo = 4.04 - 10~* Pa-s 
and 4.88 - 10~* Pa-s for theoretical harmonic and anharmonic 
c„, respectively, a = 1.22 - 10""* K~*. 



law, r/ — r]Q exp 



The data were subsequently 



T-To^ 

extrapolated on the temperature range of experimental 
Cy because the experimental was measured in a wider 
temperature range than 77, and used rj in Eqs. (jl3ll4p 
to calculate c„. We note that in Eqs. (|13ll4p de- 
pends on both rj and and is therefore sensitive to 
the viscosity fit. We find that as long as the VFT fit has 
physical values of parameters (e.g. positive To); ^ reason- 
able agreement between the calculated and experimental 
Cv is found, as discussed below in more detail. 

Notably, Eqs. ([T3)) and (fT4|) have no free fitting pa- 
rameters. Indeed, a single parameter in Eq. (IT3)) . tdGoo, 
is defined by the liquid properties. Eq. (jl4p has an extra 
parameter a, which is equally governed by liquid proper- 
ties. In practice, however, the values of td and Goo are 
not known precisely for all liquids. Therefore, we have 
varied tdG 00 around its approximate experimental values 
(see below) to fit the calculated Cv to the experimental 
data. 

In Figures ([IMl), we show the experimental Cy where 
the electronic contribution was subtracted d, Q . We also 
show Cy , calculated from both harmonic and anharmonic 
theory using Eqs. (|13|14p . Overall, Figures ([2][6]) show 
reasonably good agreement between the theoretical and 
experimental c^. 

The agreement is somewhat worse at low temperature, 
however the maximal difference between the predicted 
and experimental values is comparable with the experi- 
mental uncertainty of c„ of 0.1-0.2 J/K Q. We note that 
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FIG. 4: (Color online). Experimental and theoretical c„ for 
liquid Cs. Experimental c„ is shown in square symbols. Cal- 
culated harmonic and anharmonic Cv are shown in the bottom 
and thick black curve, respectively. tdGoo ~ 1.19 - 10~* Pa-s 
and 2.07- 10""* Pa-s for theoretical harmonic and anharmonic 
Cv, respectively, a = 4.92 - 10"* K"^. 



similar cx ^ terms appear if Eq. (IHl) is used. Here, the 
second term is small at low temperature where r ^ td 
and therefore 77 » tdGoo- The last term depends on 
both 77 and jf^- is largest at low temperature, there- 
fore the last term is most sensitive to the fitted slope of 
77 at low temperature and its extrapolation, and is ex- 
pected to contribute most to the difference between the 
calculated and experimental Cy. 



6 



300 400 500 600 700 800 900 1,000 




3.5 



^ 3 



2.5 



Experiment 
Theory: anharmonic 
Theory: harmonic 




300 



400 



500 



600 700 
T(K) 



' ' I ' 
800 



■3.5 



■2.5 



900 1,000 



FIG. 5: (Color online). Experimental and theoretical for 
liquid Rb. Experimental Cv is shown in square symbols. Cal- 
culated harmonic and anharmonic c„ are shown in the bottom 
and thick black curve, respectively. tdGoo ~ 1-46 • 10~* Pa-s 
and 1.914- 10~* Pa-s for theoretical harmonic and anharmonic 
Cv, respectively, a — 4.52 ■ 10"* K~^. 
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FIG. 6: (Color online). Experimental and theoretical c„ for 
liquid Sn. Experimental c„ is shown in square symbols. Cal- 
culated harmonic and anharmonic Ci, are shown in the bottom 
and thick black curve, respectively. tbGoo ~ 4.67 - 10"* Pa-s 
and 5.67- 10"* Pa-s for theoretical harmonic and anharmonic 
Cv, respectively, a = 1.11 - 10~* K~^. 



Importantly, the theoretical curves in Figures (l2]|6]) are 
calculated using physically sensible values of parameter 
tbGcx in the (l-6)-10~^ Pa-s range (see captions in Fig- 
ures (l2]|6l)). For example, if Goo of Hg is about 5 GPa 
[23| . the used values of tdG oo in Figure ^ imply td of 
about 0.1 ps, a typical value of Debye vibrational period 
which is, furthermore, consistent with the experimental 
value :28]. Similarly, if Goo of Cs is 0.77 GPa 0, the 



used values of tdGqo in Figure ([4]) imply td of about 
0.1 ps, equally consistent with experiments [2§|. We also 
note that the used value of TmtGoo is close to that in 
other liquids [l9| . 

Similarly to tdGoo, the values of a used to calculate 
the anharmonic c„ (see captions to Figures ^l^) are 
close to experimental ones. The experimental a for Hg, 
In, Rb, Cs and Sn are 1.8 - lO"* K'^, Lll-lO""' K"!, 
3-10-4 K^i, 3-10-4 and 0.87-10-* R-^, respectively 
in reasonably good agreement with the values used 
in Figures pie]) . 

We observe that for all five liquids, harmonic Cy is sys- 
tematically below the experimental value at low temper- 
ature. This is not surprising because the maximal low- 
temperature value of harmonic Cv is 3, according to Eq. 
(fT3| . As Figures ([2][6|) show, the anharmonic contribu- 
tion, calculated using Eq. (|14p . is important in order to 
reproduce the experimental behavior. 

In summary, we find that overall, theoretical and ex- 
perimental values of c„ agree reasonably well, using phys- 
ically sensible values of parameters tdGoo and a. 

We note that above we have discussed the behavior of 
Cy. The quantity that is frequently measured in the ex- 
periment is the constant-pressure heat capacity, Cp. For 
several low-viscous liquids, Gp weakly depends on tem- 
perature, in contrast to the constant- volume heat capac- 
ity, Gv [1, [Hi [s^l - On the other hand, for highly viscous 
liquids, Gr, increases, approximately linearly, with tem- 
perature |3l| . This behavior can be rationalized in the 
proposed picture as follows. We re- write the known re- 
lationship, Cp^ Gy + VTa^B [H, as Cp = c„(l laT), 
where Cp = ^ and 7 = is the Griineisen pa- 

rameter, which is on the order of unity in various sys- 
tems [32|. For highly viscous liquids where t 3> td, 
the last term in Eq. (fTTj) representing the contribu- 
tion from the decreasing number of transverse waves, 

3 - (^)^^ > can be neglected. Then, c„ = 3(1 + aT), 

giving Cp = 3(1 + (7 + l)aT), where we have neglected 
the square of the small parameter aT . This is consistent 
with the linear increase of experimental Cp 3lJ. On the 
other hand, in low- viscous liquids where r approaches 
Td at high temperature, the contribution of the decreas- 
ing number of transverse modes in Eq. (|11|) can not be 
neglected. In this case, Cp — Ci,(l -I- ^aT) is the prod- 
uct of temperature-decreasing Cy due to the decreasing 
contribution of transverse waves (see Figures ([2][6])) and 
temperature-increasing term, oc aT. As a result, Cp is 
less sensitive to temperature, as is seen experimentally. 



TWO APPROACHES TO LIQUID ENERGY: 
DISCUSSION 

It is interesting to compare the proposed approach of 
calculating liquid energy to previous theories. A tradi- 
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tional way to calculate liquid energy was based on the 
calculation of liquid potential energy in addition to gas 
kinetic energy In this sense, the approach is from 

the gas state, in that a liquid is considered to be a gas 
with interactions switched on. In the approach from the 
gas state, Ikiuid energy can be written in a generalized 
form as lll- 



E = ^NT + I UFdV 



(15) 



where U is the interaction energy between atoms and F 
is the appropriately normalized correlation function. 

In the proposed approach from the solid state, liquid 
energy is calculated on the basis of solid-like phonons, 
which therefore accounts for strong interactions from the 
outset. For convenience and reference to the foregoing 
discussion, we re- write the harmonic energy, Eq. ([7]), as 
a function of r: 



(16) 



Early approaches based on Eq. (jT5|) were most success- 
fully used for dilute or weakly-interacting systems 
For example. Van Der Waals equation serves to illustrate 
the effects of interactions and excluded volume. At the 
same time, c„ of the Van Der Waals liquid is | , the same 
as that of the ideal gas [H, stressing that this system is 
far from a real dense liquid. 

Historically, the same general approach from the gas 
state was also used to discuss real dense and strongly- 
interacting liquids [H-H, 13, 14, 33 1. We now propose 
that the approach to liquids from the solid state offers a 
number of important advantages. 

First, experimental liquid Cy changes from 3 right af- 
ter the melting point to about 2 at high temperature (see 
Figures ([2][6])). In the approach from the gas state, the 
second term in Eq. ([TSl) should be as large as the first 



one {^NT) in order to reproduce the low-temperature 
experimental behavior. On the other hand, the starting 
value of Ci, = 3 in Eq. (jl6[) is already equal to the exper- 
imental one, suggesting that this equation offers a better 
starting point to discuss liquid c„. 

Second, the approach from the gas state does not read- 
ily explain the experimental decrease of Cy from 3 to 2. 
The increase of temperature usually results in broadening 
of correlation functions, but without the specific form of 
F at each temperature as well as U and explicit integra- 
tion, it is not easy to see how Eq. (fT5|) explains the exper- 
imental behavior. On the other hand, Eq. ()16|) explains 
the experimental data of Cy in a simple and straightfor- 
ward way. 

Third, both harmonic and anharmonic effects in the 
approach from the gas state are contained in U and F, 



and their separation is not straightforward in general. On 
the other hand, the anharmonic effects are readily calcu- 
lated in the approach from the solid state in addition to 
harmonic terms, from the change of phonon frequencies 
and thermal expansion (see Eq. ([TTt '). Consequently, the 
experimental increase of Cy above the Dulong-Petit value 
at low temperature is attributed to anharmonicity. 

Fourth, it is important to stress that F and U are not 
generally available for liquids, apart from simple systems 
such as hard spheres, Lennard- Jones liquids that model 
Ar, Kr and so on. For these systems, F and U can be 
determined from experiments or simulations and subse- 
quently used in Eq. (fTS]) to calculate the liquid energy. 
Unfortunately, neither F nor U are available for liquids 
with any larger degree of complexity of structure or in- 
teractions. For example, many-body correlations fl3, 14 1 
and network effects can be strong in familiar liquid sys- 
tems such as olive oil, Si02, Se, glycerol or even water 
[i^ l , resulting in complicated structural correlation func- 
tions that can not be reduced to simple two- or even 
three-body correlations that are often used in Eq. ([T5|). 
As discussed in Ref. Q, approximations become diffi- 
cult to control when the order of correlation functions 
already exceeds three-body correlations. Similarly, it is 
challenging to extract multiple correlation functions from 
the experiment. The same problems exist for interatomic 
interactions which can be equally multi-body and com- 
plex, and consequently not amenable to determination in 
experiments or simulations. 

On the other hand, r is available much more widely, 
and is readily measured in experiments of various types, 
including dielectric relaxation, NMR and so on. Equally, 
viscosity rj = GooT is routinely measured in liquids using 
well-known methods, no matter how complicated struc- 
tural correlations functions or interatomic interactions in 
a liquid are. Therefore, Eq. (|T6l) can be more readily used 
to calculate liquid energy from the practical perspective. 
In future work, we plan to calculate heat capacities of 
complex liquids in the proposed approach. 

Importantly, Eqs. (ITSI) and are related, in that r 
that enters Eq. ([T6)) is itself governed by F and U in Eq. 
([T5|) : T — t{F,U). However, this function is generally 
very complex because it is defined by the activation bar- 
riers set by F and U and not by their equilibrium proper- 
ties. Calculating r for a given set of F and U can not be 
done analytically using a general procedure because the 
potential energy landscape gets complex even for a small 
number of particles and possible configurations. This 
task can be done in a simulation but is equally challeng- 
ing for realistic systems because F and U are not known 
except for a small set of simple systems. At the same 
time, going to the deeper level and calculating the en- 
ergy from first-principles using F and U is not required 
in our approach because it is r that defines the liquid 
phonon states from the physical point of view, as pro- 
posed by Frenkel originally 2]. Consequently, F and U 
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do not feature in the approach from the sohd state (see 
Eq. (IT6)) '). just as they do not in the energy of a sohd. 

Fifth and finahy, the proposed approach from the sohd 
state can be said to be more general, or universal, as com- 
pared to the approach from the gas state, in the follow- 
ing sense. In the approach from the gas state, the energy 
strongly depends on U as well as on F (see Eq. p3|) . It is 
for this reason that Landau and Lifshitz state the liquid 
energy strongly depends on the specific form of U and 
therefore can not be calculated in general form [H. Lets 
now consider liquids with very different structural corre- 
lations and interatomic interactions such as, for example, 
H2O, Hg, AsS, olive oil and glycerol. Even if their energy 
could be calculated on the basis of Eq. ([T5|) , large differ- 
ences in U and F would imply different values of energy, 
in agreement with the argument of Landau and Lifshitz 
On the other hand, as long as r of the above liquids 
is the same at certain temperature (different for each liq- 
uid), Eq. (fT6|) predicts that their energy is the same. In 
this sense, expressing the liquid energy as a function of r 
only is a more general description because r is common 
to all liquids whereas F and U are system-specific as in 
the above examples. 



CONCLUSIONS 

In summary, we have developed the approach to liq- 
uids from the solid state by incorporating the effects of 
anharmonicity and thermal expansion. In this approach, 
we have found a good agreement between theoretical pre- 
dictions and experimental data for five commonly dis- 
cussed liquids. We have proposed that the approach to 
liquids from the solid state offers a number of advantages 
as compared to the previous approach from the gas state. 

We are grateful to V. V. Brazhkin for discussions and 
to SEPnet, EPSRC, Myerscough Bequest and the School 
of Physics in Queen Mary University for support. 
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